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Abstract. The first measurement of is presented. At low x the measurement 

■ indicates a non-zero hadron-like component to i^ c . At large x the measurement con- 
stitutes a test of perturbative QCD at next-to-leading order, with only m c and a s as 
free parameters, with a precision of 0(40%). 
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O : INTRODUCTION 

(N 

For about 20 years measurements of photon structure functions give deep insight 
into the rich structure of a fundamental gauge boson, the photon. A recent review 
on this subject can be found in [1]. Here, the discussion is restricted to the mea- 
surement of F^c which recently has been achieved for the first time. Only the main 

■ features of the analysis are given, the experimental details can be found in [2]. 
The differential cross-section for deep inelastic electron-photon scattering, shown 

<D in Figure 1, is given by 
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(l + (1 - yf) F2(x, Q 2 ) - y 2 F?(x, Q 2 )] . (1) 
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b : 

Here Q 2 is the absolute value of the four momentum squared of the exchanged vir- 
tual photon, 7*, x and y are the usual dimensionless variables of deep inelastic scat- 
tering and a is the fine structure constant. In experimental analyses y 2 is usually 
small. Consequently, the term proportional to the longitudinal structure function 
Fl can be neglected and the differential cross-section is directly proportional to 
F2 , which is related to the sum over the quark parton distribution functions q 1 of 
the quasi-real photon, 7, via 

Fj(x, Q 2 ) = x |f < [q\x, Q 2 ) + if(x, Q 2 )} . (2) 

q=u,d,s 



^ Invited talk given at the PHOTON 2000 Conference, Ambleside, UK, August 26-31, 2000, to 
appear in the Proceedings. 
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FIGURE 1. Examples of leading order diagrams contributing to (left) the point-like, and (right) 
the hadron-like part of F^c- 



Due to the large scale established by their masses, the contribution to F2 from 
heavy quarks can be calculated in perturbative QCD. At present collider energies 
only the contribution of charm quarks is important. Like the structure function 
for light quarks, F% c receives contributions from the point-like and the hadron-like 
component of the photon shown in Figure 1. 

Because of the charge of the charm quarks their contribution to F2 is large 
and the importance increases for increasing values of Q 2 , as can be seen from 
Figure 2, which shows the contributions from light quarks and from charm quarks 
separately, as predicted by the GRV parametrisations [3]. Charm quarks can only 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X 

FIGURE 2. The structure function F% for u,d,s quarks alone and for u,d,s,c quarks, as a 
function of x and for two different values of Q 2 . 

be produced if the photon-photon invariant mass W is at least twice the mass of the 
charm quarks m c . Using x = Q 2 /(Q 2 + W 2 ) this leads to the varying production 
threshold in x as a function of Q 2 seen in Figure 2. 

Close to the production threshold, the point-like contribution to is accu- 
rately approximated by the prediction of the lowest order Bethe-Heitler formula. 
For quasi-real photons also the next-to-leading order (NLO) predictions have been 
calculated in [4] . For the hadron-like contribution the photon-quark coupling must 
be replaced by the gluon-quark coupling, and the Bethe-Heitler formula has to be 
integrated over the allowed range in fractional momentum of the gluon using a 
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FIGURE 3. The predictions of the SaSID (full) and the GRS (dash) parametrisations for 
Q 2 = 30 GeV 2 and for (a) P 2 = and (b) P 2 = 1 GeV 2 . In addition, the point-like contribution 
for m c — 1.5 GeV(dot-dash) is shown. 



parametrisation of the gluon distribution function of the photon, see e.g. [1]. 

The predicted behaviour of the point-like and hadron-like component of K7 C 
for different values of m c is shown in Figure 3, using the SaS ID [5] and GRS [6] 
parametrisations. For x > 0.1 the structure function is saturated by the point-like 
component which is only slowly suppressed for increasing virtualities P 2 of the 
quasi- real photon. In contrast, the hadron-like contribution dominates at small 
values of x and decreases much faster for increasing P 2 . Finally, lowering the mass 
of the charm quarks leads to a higher threshold in x. 

Given this predicted behaviour, the region of x > 0.1 can be used to test a purely 
perturbative NLO QCD prediction with only m c and a s as free parameters, and 
the low x behaviour mainly probes the gluon distribution function of the photon. 

MONTE CARLO MODELS 

The LO Monte Carlo generators HERWIG 5.9 [7] and Vermaseren [8] are used, 
both with m c = 1.5 GeV. In HERWIG charm production is modelled using matrix 
elements for massless charm quarks, together with the GRV parametrisation for 
the parton distributions of the photon, again for massless charm quarks. The 
effect of the charm quark mass is only accounted for by not simulating events 
with W < 2m c . Due to the massless approach used in HERWIG and the crude 
treatment at threshold, the predicted charm production cross-section is likely to be 
too large. The Vermaseren generator is based on the Quark Parton Model (QPM) 



and consequently does not take into account the hadron-like component of the 
photon structure. However, the complete dependence of the cross-section on the 
different photon helicities is modelled. 

THE MEASUREMENT OF F% c 

The measurement of proceeds along the same lines as the measurement of K7 
with the addition of the identification of the charm quarks via the reconstruction 
of D* mesons. 

Events are selected with an energy of the scattered electron above 50 GeV, mea- 
sured in the angular ranges (a) 33 — 55 mrad or (b) 60 — 120 mrad from either 
beam direction, thereby covering the approximate range in Q 2 of 5 — 100 GeV 2 . 
The visible hadronic mass W vis is required to be below 60 GeV. Charm quarks 
are identified via D* — > D°7r, followed by D° — > Ktt or D° — > Ktttttt. Using 
f(c — > D*) = 0.235 ± 0.011 and the branching ratios of the D° decay modes of 
0.02630 ± 0.00082 and 0.0519 ± 0.0029, this analysis covers only about 4% of all 
events containing a pair of charm quarks. For a clear acceptance the D* mesons 
are further required to fulfill, \i] D *\ < 1.5 and p^* > 1 or 3 GeV for (a) or (b). 
Together with a typical selection efficiency of about 25% only about 1% of all cc 
events are positively identified. 
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FIGURE 4. Mass difference A M for the anti-tagged and tagged sample. 

Figure 4 shows the distribution of the difference between the D* and the D° 
candidate mass. In both samples a clear peak is visible around AM = 145.4 MeV. 
Subtracting the background, obtained from a fit to the upper sideband of the signal, 
29.8 ± 5.9(stat) D* mesons are found in the peak region. 

Figure 5 shows the distributions of W V - 1S and of the measured Q 2 in comparison to 
the predictions of the HERWIG and Vermaseren Monte Carlo generators normalised 
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FIGURE 5. The distributions of a) the visible invariant mass W v i s , and b) the negative 
four- momentum squared Q 2 . 



to the number of data events. Both Monte Carlo generators give a good description 
of the shape of the data distributions. 

The cross-section for D* production is determined in the well-measured kine- 
matic range described above. Based on this and the extrapolation factors obtained 
from the Monte Carlo models the full cross-section <r(e + e~ — > e + e~cc) and F£ c are 
evaluated in two bins of x with 0.0014 < x < 0.1 and 0.1 < x < 0.87. 

For x > 0.1 the predictions of the Vermaseren and HERWIG Monte Carlo models 
are very similar. In contrast, for x < 0.1, there are two main differences. Firstly the 
selection efficiency for events with D* mesons fulfilling the kinematical requirements 
is different, and secondly, and even more important, the predicted cross-section 
within the invisible phase-space is largely different for the two models, resulting in 
different extrapolation factors, 12.9/5.1 for HERWIG/ Vermaseren. 

Since the hadron-like contribution is neglected in the QPM, the Vermaseren cross- 
section is much smaller than the LO and the NLO cross-section for x < 0.1. In 
contrast, mainly due to the massless approach taken, the prediction from HERWIG 
is higher than the cross-section from the LO and the NLO calculation. Therefore 
it is likely that the correct cross-section, and therefore the correct extrapolation 
factor, lies within the range of the two Monte Carlo predictions. 

The measured K7 C for (Q 2 ) = 20 GeV 2 is shown in Figure 6. The central values 
are obtained by averaging the results using the HERWIG and Vermaseren Monte 
Carlo models, and half the difference is taken as extrapolation error, which domi- 
nates the uncertainty for x < 0.1. The NLO prediction is based on m c = 1.5 GeV, 
the renormalisation and factorisation scales are chosen to be /xr = /xf = Q, and the 
hadron-like contribution to F^ c uses the GRV parametrisation. The NLO correc- 
tions are small for the whole x range. The band for the NLO calculation is evaluated 
by varying m c between 1.3 and 1.7 GeV and using Q/2 < /xr = /xf < 2 Q. 

For x > 0.1 the error of the measured cross-section is dominated by the statistical 
uncertainty, and the NLO calculation with only m c and a s as free parameters is in 
good agreement with the data. In contrast, for x < 0.1, the result suffers from the 
strong model dependence discussed above. Despite this uncertainty the corrected 
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FIGURE 6. F^ c compared to several predictions explained in the text. 



data suggest a cross-section which is above the purely point-like component, i.e. 
the hadron-like component of K7 C is non-zero. 

Conclusion and Outlook 

In conclusion, for x > 0.1, the purely perturbative NLO calculation is in good 
agreement with the measurement and for x < 0.1, the measurement suggests a 
non-zero hadron-like component of F^ c . 

By using the massive matrix elements available in HERWIG6.1 and the full inte- 
grated luminosity of more than 500 pb" 1 of the LEP2 programme, the measurement 
is likely to be improved considerably, both concerning the statistical and the sys- 
tematic error. 
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